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1. INTRODUCTION 

The maximum power generation of a wind turbine has been interested in several decades and many algo- 
rithms have been suggested. According to [1], previous maximum power point tracking (MPPT) methods can be 
listed into three groups including indirect power controller, direct power controller, and others. The indirect power 
controller which aims to maximize mechanical power by using tip-speed ratio [2], [3], optimal torque [4], and power 
signal feedback MPPT algorithm [5] is simplicity and only allows the wind turbine to track its MPPT-curve quickly 
when a wind speed measurement is precise and instantaneous. In the case of an unavailable wind measurement, a 
wind turbine using the indirect power controller fails to track its maximum power point quickly and accurately [6], 
[7]. The direct power controller which maximizes electric power by using the perturbation and observation (P&O) 
algorithm such as: Hill climb search [8], incremental conductance [9], optimal-relation based MPPT algorithm [10], 
[11], and hybrid MPPT algorithm [12] does not require any wind turbine knowledge and available anemometer. Un- 
fortunately, a wind turbine using the direct power controller cannot track its maximum power because this controller 
cannot recognize instantaneously the variation in wind speed. Until now, this controller has been implemented to 
adjust the voltage and current of DC circuit in a permanent magnetic synchronous generator wind turbine. The last 
group is developed based on soft computing techniques like Fuzzy [13] and Neural network [14]. A wind turbine us- 
ing these methods only has a good performance when the full information of the wind turbine is available. However, 
these methods are complexity and large memory requirement. Hence, a new MPPT scheme for DFIG wind turbines 
should be researched. 

To control the generator-wind turbine, proportional-integral (PI) control is normally implemented because 
of its simplicity [7], [15], [16]. However, by using the PI control, we cannot guarantee the wind turbine system will 
become stable operation [17], [18]. Recently, control laws based on sliding mode were suggested for rotor speed 
adjustment [19], [20]. However, these control laws require an available wind speed measurement. Hence, we need 
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Fig. 1. DFIG wind turbine configuration 
to propose a new control law for rotor speed in the DFIG-wind turbine. 

In this research, we suggest a new scheme to extract the maximum available power from a wind turbine 
employing DFIG. This scheme is developed from the feedback power algorithm but in this research, we do not 
require an anemometer. New control laws which are developed upon Lyapunov function for rotor speed, current, and 
voltage are proposed. This scheme is validated through numerical simulations of a wind turbine employing DFIG. 
From simulation results, we will analyze and compare with the simulation results of a wind turbine using an old 
MPPT scheme. 


2. DFIG WIND TURBINE 
A DFIG-wind turbine is described in previous publications, as shown in Fig. 1. Generally, it consists of a 
wind turbine, shaft-gearbox, doubly-fed induction generator (DFIG), and back-to-back converter. 


2.1. Wind turbine 


When the wind turbine is rotating at a speed of w, and wind speed at the wind turbine is V,,,, the mechanical 
power of the turbine is calculated through blade length R, air density p, and power coefficient C,(A, 6) [21] 


| 
Pm(t) © spr R°Op(A, B)Vo(t)- (1) 


The wind turbine’s power coefficient C,,(A, 3) depends on both pitch angle £ and tip speed ratio [16] 


2 Rw, (t) 


sg 


(2) 


At a constant A, when £ increases C’,(/) will be decreased. In contrary, at a constant 3, C'p(A) reaches to a maximum 
value Cp(Aopt) at ÀA = Aopt- 


2.2. DFIG 


The DFIG’s main objective is to convert the mechanical power P,,, on the wind turbine shaft to electricity 
power P. Relationship between P, and Pm is described through the DFIG-wind turbine’s inertia J 


Jur(t) Sarl) = P(t) — Pelt) 6) 


Generally, the DFIG is an induction generator so its rotor slip is defined as 


a(t) E AN (4) 


Ws 


where pn is the number of pole pairs; N is the gearbox ratio; ws is the rotational speed of stator flux. 
T T 
In dq frame, the stator voltage v, £ [usa Usg] and rotor voltage v, £ [vra Urq] are computed from 
i ; a BL : ; * Sale T 
the stator current i, 2 lisa isq] , rotor current i, £ lira irq] , stator flux Y, (t) = [Wsa(t) Veq(t)] , rotor 


T . ; : : 
flux Y, (t) = [Wra(t) Urg (t)] , rotor resistance 7,, rotor inductance Ly, stator resistance rs, stator inductance 
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L., magnetizing inductance L, as [22] 


d 
+w,OW,(t) + a Vslt) 


v,(t) = rri (t) + s(t)wsO (P, (t) + Lao, (5) 
W(t) = Lsis(t) + Lmir(t) 
W,.(t) = Lrir(t) + Lmis(t) 


aJo =l 
where © a |} Ha 


Lemma 1. If we neglect the stator resistance, rs = 0, and choose the dq frame so that Y(t) = [Wsa 0] E then we 
can write the state-space equation of the DFIG (5) as 


Li = A(O) + ov) + dlt), © 
where 
A Te A = A Lm 
Tl- F Ali) -or —wrs(), d(t) &-*s(0) [o] IDE k i (7) 


Proof. Obviously, if we ensure Y(t) = [Tsa 0] T then by using (5), we have 


W(t) = Lsig(t) + Lmip(t) = [Ysa 0] |, “ w(t) =06 Li) = —Lim in( 8) (8) 
By using (8) and rs = 0 in (5), we have 
v.(t) =w.OW,(t) = [0 wss] =[0 V], (9) 
where we used V; = ||vs(t)|| = |wsWsa|. From (8) and (9), we have 
it) =- Fi) + |G]. POPO (10) 
From (10) and (5), we have 
W(t) =i) + P* [Ve 0], S(t) = o Tilt), a1) 
where we used o = Ly — EA J Ls. To use (11) and (5), we have 
elt) =rrin (t) + owes (OCE) Hw OEE [Ve 0)” +oSi,(0, 2) 


Svs 


From (12), we can extract (6) easily. 


From [22] and by using (9) and (10), we calculate the stator side active power P, in the DFIG as 


Lim, . 
P,(t) = Usaisa + Usqtsq = — 5, Veir). (13) 


2.3. Grid side converter 


For the DFIG, the electricity frequency on the rotor side always depends on the rotor speed wp. Hence, 
to interface to the connected grid, a back to back converter which includes a rotor side converter(RSC), a DC-link, 
and a grid side converter (GSC) must be installed on the rotor side of the DFIG [23]. Normally, to reduce harmonic 
components generated by the GSC, a filter which consists of a resistor Ry, an inductor Ly in series and a power 
factor correction pf in parallel is used as Fig. 1. 
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In dq frame, the relationship of voltage v, = [vga vga] and current ig = liga iga] of the GSC is 
written [24] 


d, . 1 
qr = Agig(t) + dgy + te (14) 
where 
A, =—L;'Rylz—w,@, dg =—L;1[Vs 0]. (15) 


3. MAXIMUM WIND POWER EXTRACTION SCHEME 
The optimal power control region of a wind turbine is limited by [25] 


D = {(w,, Vw) | Wrmin < Wr < Wrrated; Vwmin < Vw < Vuwrateds B = 0, and Cp(à, p) > 0}, 


where Wrmin and wrrateq are the minimum and rated rotor speed, respectively; Viymin and Vwratea stand for the 
minimum and rated wind speed; ( is the blade system’s pitch angle. Hence, when the wind turbine operates in D, 
the tip-speed ratio and the rotor speed reference are limited by 


A RwWrmin 
Amin a) a 


< X(t) <Amax + max{r| Cp(A, 8) > 0}, 
Vwrated 


Wrmin <Wrref < Wrrated- 


From (1), to extract the maximization of the mechanical power, we must adjust w, to obtain the maximiza- 
tion of C,(A(w,, Vw)). For any wind turbine, we have [24] 


Cmax = Cp(Aopt), Aopt = arg max Cp(A). (16) 
As the wind turbine operates at Apt, the rotor speed becomes the optimal rotor speed 
Xo Vw 
Wropt (Vw) = a (17) 
and the mechanical power becomes maximal 
1 
EA Pm (wr, Vw) z 5 PTR Cpmax Vi = koptWropt (Vw); (18) 
1 C 
kont = -pr R? a, (19) 
ae” Nopi 
In this paper, the power coefficient in D is given as 
165.2842 
Cp(à) = (= — 16.8693) ek + 0.009A, (20) 


it has an unique maximum point of Cymax = 0.4 at Aopt = 6.7562, and its mechanical power at different wind speed 
as shown in Fig. 2. 
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Fig. 2. MPPT curve of wind turbine 


Remark 1. From (1), (2), and (19), we have [24] 
Pm (t) = kopiis (t) = C(t), (t) (Wropt (t) — Wr (t)), (21) 


pelt) Se MO ~ C0) 
2X(t) A(t) — Acpt 


C(Wr, Vw) = Tp > 0. (22) 
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3.1. Conventional MPPT-curve scheme 


The conventional MPPT-curve method [7] makes 


Wr > Wrref = \/ P./Kopt- (23) 


The problem is this conventional method cannot track quickly the maximum power point. Hence, we need to propose 
a new method as next subsection. 


3.2. Proposal of maximum available power scheme 


The subsection aims to propose a new scheme to maximize P,,,(w,, Vw) or minimize the error |Wropt — 
wr(t)|. To obtain this target, we propose wres (t) satisfying 


d d 
koptWrrep (t) a Pe (t ) + akopt (w Were f (t) TW z(t )) A hwy (t +e (9Ymax — hwy (t )) (24) 


where k, a, Ymax are positive constants; and 


0 for |k— <2) < ymax 
q= 1 for kw; (t) > Ymax (25) 


-1 f k max: 
or gÀ < =y 


From (23) and (24), we can see that the proposed scheme is developed from the conventional MPPT method. 


4. CONTROLLER DESIGN FOR DFIG 
4.1. Rotor side control 


The purpose of RSC controller is to reduce the errors (ira — irdref) and (wr — wrref ) in which iraref and 
Wrref are the reference of irq and w,, respectively. From (3), (6), (13), to make w, converge tO wrref, we can adjust 
irq to a reference value trqref corresponding wrref. In [7], this task is carried by traditional PI controls. In this 
research, to obtain the above target, we design v, of the DFIG (5) as 


vr(t) = A (t)i (t) + d(t) — oK, (irre (t) — ir (t)) — oT isali, (26) 
irref (t) = [irdret (t) irq (t) + kprewrret (t) + kir j: Cwrref (r)dr] 7 ) (27) 
Cvdrrat (E) = Wires (t) — Wr (t), (28) 


where, kpr > 0, kir > 0 and matrix K, > 0. 


Theorem 1. When the DFIG-wind turbine operates in D, if wrref and v, of the DFIG (5) are designed as (24) and 
(26), respectively and if there exist positive constants u1, 42, and b satisfying 


min (xe ja EINA s) >bmð, (29) 
Hı H2 

min (2k - pi2€(t) — yio 1] K,K] [0 1) > buikpr, (30) 

K, +K] - yK] [0 1]' [0 1] K, > bm, (31) 
where 

FAI—2w1-@)), E e — FZ, e@) 4 (1 - a) A 
xy di Pt w(t)’ 0 

then 


maxw, |y) 


lim (irref (t) — ir(t)) = 0, jim (weet (t) = w(t) =0, |Wropt — wr| < Tb 
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Proof. Let define 


e(t) = [Crmar(t) eir(t)]" = [w3ee(t) —w3(t) irres(t) — i, (2)] | 


em(t) = [europe lt) €t] = [Wope(t) —wr(t) wi ee(t) — w(t) irre g(t) — i, (8)] 


By substituting (26) into (6), we have 


dz ; : ; 
dt (ipret (t) — irp (t)) = —K, (ipret (t) =li, (t)) . (32) 
(27) can be rewritten as 
d d 
lyr erat) = kireva (t) + [0 1] Feil), (33) 
and by substituting (32) into (34), we have 
ee (t) = —kirewnes (t) — [0 1] K,eir(t), 84) 
Then, 
EL e, (t) = -Q,e,(t) (35) 
rdt T a i hae of $ 
_ |kpr 0 lR [0 1] K, 
E, = L 1 > 0, Qr = kK K, bi (36) 


When we define a Lyapunov function as V, = e (t)E,-e,(t), its derivative is 


ART d z 
a =e, (Er ser (t) + (50) E,e,(t). (37) 


By substituting (35) into (37), and noting that Q, + Q, = Q,, we have 
d 
dt 

Furthermore, by substituting (24) into (3) 


V, = —e} (t) (Qr + Q; ) e,(t) = -e! ()Q,e,(t) < -Amin (Qr)e; (t)e,(t). (38) 


d 


Jur (t) Šilt) = Prat) — kopttirotlt) + kop (Seat) — W) + KL t2(t) + a? (auma — BSu) 
= Pra (t) — hope) + koplo? (t) — fros (t)) + ket) + akoni erlt) — 02(0) 
+ q? (Ymax — bu? ) 
= COJE erg t) + E(B or are lE) + Umax + KC Ber (ter), G9) 


where we use (21) and (a — 1)Kopt (w3 p(t) — w3 (t)) = E(t)wy (t) (Wrret(t) — w,(t)). By using J & J —2k(1—q?), 
(39) becomes 


ad = Pas 

I Slt) = euron 8) + Eeun lt) + ase (40) 
It means 

~d 

J gorom E) = =C) wrons (t) = Elb) ewr lt) + V(t), (41) 

> d CYmaz 
where we used y(t) = I(T Wropt o ). Hence, 
d 
Em gen) = Qm (t)em (t) + Mm (t) (42) 
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Saa CO EO Onz 
Emn = l ral , Qm (t) = 0 kir [0 1] K,| , 
3x1 í 0271 0271 K, 
We define Lyapunov function as Vm = €m (t) | Em€m/(t), its derivative is 
d eea ad d z 
qe = e,,, (t) En 7 em(t) + (Enge) eEm(t). (43) 
By using (42) in (43), we have 
d 
= —em(t)' Qm (t)er (t) + em(t) Mm (t) — em(t)' Qm (t)! €m (t) +Mm(t)' e€m(t). (44) 
Noted that for pı > 0 
Em (t) Mm (t) + M(t) "@m(t) = 2ewropi (IE) < Coron OIEI + ealy@)| 
< em (t) Mmi (tem (t) + mlr) (45) 
where we used Mmı(t) = diag(|y(t)|/ u1, 0,0, 0). Hence, 
d - 
—Vm S —em(t)' Qm(tem(t) + paly(t), (46) 


dt 


where Qm = Qm (t) + Qm(t)" — Mmı (t). 
Remark 2. For (29)-(31), with u2 > 0, we have 


Hı 


EN si (2 KOl EO a me ie 1] K,KT fT] Ke + Ky -je fil (0 i1x,) 


and certainly, Q, (t) > 0. 


Hence, according to the Lyapunov Stability Theory, (38) and (46) give us limy_,.0. €r (t) = 0 and |wropt — 


wel <Vn/F < ymax. WOI. 


4.2. Grid-Side Control 


In this section, we propose a new control law such that Vgc and 74, are maintained at their references Vacret 
and tggref, respectively. To maintain Vg, at Vacrer, we need to make i,q converge to igaref Corresponding to Vacret. 


Theorem 2. For any Vacref and igqref, if Vg of the GSC (14) are designed as 
d., g T 
v(t) =L; pier (t) + K,ei,(t) — Agi,(t)} + [Vs 0] 
r ; . T 
lgr (t) = liga(t) + kpg€v (t) ale kig J Ey (r)dr igarett)] 
€v (t) >. Veret (t) ii V(t), Cig (t) = igr(t) pi ig (t) 


and if there exist kpg > 0, kig > 0, and K, > 0 with 


2 | 2kig [1 AK] , 


Q%= kī [ji o0" KT4+K, 


then 
Jim (Vierer (t) z vVi-(t)) = 0, jim (igaret (t) D tgq (t)) = 0. 


(47) 


(48) 
(49) 


(50) 
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Proof. Let define e,(t) = |e,(t) ejg(t $ . By using (47) and (14), we have 
f g g y g 


Telt) =—K,e;,(t). (51) 
Furthermore, by taking time derivative of (48) and then using (51), we have 
bya $e = —kigev(t)+ [1 0] Leili) = —kigev(t)— [1 0] Kyeig(t). (52) 
Hence, 
B, “e,(t) = —-Q,e,(t), (53) 
where 
z= [hy 8] >o a= fie D E, es 
If we use a Lyapunov function as V} = e,(t) ' E,e,(t), its time derivative will be 
d d 5 
T = e,(t)' Ey Gea!) (£, Se) e,(t) (55) 
By substituting (53) into (55), we have 
T a —e,(t) ' Qgeg(t) a e,(t)' Q; ex(t) = —e,(t) ' Qye,(t) < —Amin(Qg)eg(t) ' eg (t). (56) 


d 
Hence, if (50) holds, then a9 < 0 for all nonzero e, . This completes the proof of Theorem 2. 


5. SIMULATION RESULT AND DISCUSSION 


To evaluate the performance of the suggested MPPT scheme, we compare the simulation results of the 1.5 
MW DFIG wind turbine with that using the conventional MPPT-curve scheme with traditional PI controls [7]. In 


this research, the generator and turbine parameters [22] as shown in Table | are used. 


Table 1. Parameters of wind turbine and DFIG[22] 


Name Symbol Value 
Rated power P 1.5 MW 
The length of blade R 35.25 m 
Rated/minimum rotor speed Wprated /Wrmin 22/11 rpm 
Rated wind speed Vaated 12 m/s 
Rated stator voltage V; 690 V 
Rated stator frequency f 50 Hz 
Number of pole pairs Pn 2 p.u 
Rotor winding resistance Te 2.63 mQ 
Stator winding inductance Ls 5.6438 mH 
Rotor winding inductance Lr 5.6068 mH 
Magnetizing inductance Lim, 5.4749 mH 
Inertia of system J 445 ton.m? 


With the power coefficient (20), the region D is 


1.15 < wp < 2.3, 1.15 < wmef < 2.3, 5< Vy <12, 3.4=Amin <A < 10.239. 


In this region, (wr, Vw) as Fig. 3, which gives the minimum value ¢, min (wr, Vw) = 1.271 x 10°. 
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Fig. 4. Simulation results: ( a) wind profile; ( b) error between wropt and wp; ( c) power coefficient Cp; ( d) error 
between Pmax and Pm; (e) electrical energy output 


10 
3 Vive (m's) 


Here, we use RSC controller’s parameter as kir = 0.4J, kpr = 0.65J, K, = J diag(0.5,1), k = 0.3J, 
a= 0.2, Ymax = 0. 1Wrrated- 


For a wind profile with |2 Vo] < 0.44m/s?, the boundary of |wropt — wr| is determined as Table 2. 


Table 2. Limitation of |wropt — wr| as fi = 1 and u2 = 2 


Object q=1 q=0 
J 4.45x10° 3.12 10° 
max ¿(t) 0.8673x10° 0.8673 10° 
max |¥(t)| 0.263 x 10° 0.263 
b = 0.5915 0.1895 0.5915 
|Wropt — Wr|  1.2248rad/s 0.3775 rad/s 


When the wind profile as Fig. 4a is used, simulation results are demonstrated in Fig. 4. Fig. 4b shows that 


when the wind speed has an insignificant change, the turbine speed is almost kept up at its optimal value. Since the 
wind turbine’s large inertia, the turbine speed fails to respond instantaneously to the rapid change of the wind; this 
makes the turbine speed impossible to keep up at its optimal value. Therefore, the error |w,opt — wr| increases when 
the wind velocity changes rapidly. However, comparing with the turbine using the old MPPT scheme, by using the 
suggested scheme, the turbine speed can retain its optimal value more promptly because of the decrease in inertia 
from J to (J — a) and the error |wropt — w,| is smaller. As a result, during a rapid change in wind conditions, in the 
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Fig. 5. Error between reference signal and actual output in the controller: ( a) errors of RSC controller and ( b) errors 
of GSC controller 
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Fig. 6. Simulation results with 5% measurement noise : (a) power coefficient; ( b) error between Pmax and Pm; (©) 
electrical energy output 


proposed method, Cp restores Cymax more quickly, as shown in Fig. 4c. Obviously, by implementing the old MPPT 
scheme, the C, can be reduced to to 0.363 while by implementing the proposed scheme, this data is 0.393. 


Figure 4d shows the efficiency of the suggested method comparing with the conventional one in terms of 
mechanical power. When the wind velocity varies insignificantly, the error (Pmax — Pm) in the wind turbine using 
the suggested scheme is like that using the conventional one. However, this error becomes significant when the wind 
changes suddenly; by using the new method this error is significant smaller comparing with that using the old one 
thanks to the restoration of Cp. 


Fig. 4e indicates that to increase the turbine velocity in the period of 20s-40s, the turbine using the offered 
scheme requires a higher mechanical power comparing with that using the old one. However, the stored mechanical 
power is returned in the interval of 60s-75s in which the rotor speed decreases. Hence, in the period of 20s-40s, 
comparing with the DFIG using the suggested scheme, the electric energy generated by the DFIG using the old 
MPPT method is little higher but in the 60s-75s interval, it becomes opposite. As a result, accumulating to the end 
of simulation, the wind turbine using the conventional method fails to generate the electrical energy in total as high 
as that using the proposed method, as Fig. 4e. This indicates the quality of the suggested MPPT scheme. 


Fig. 5 shows the control quality of the RSC and GSC. Both (w3 ef — w3) and (iraref — ira) in Fig. 5a are 
very small, it means w, and tyrą track their reference values; in other words, the control law proposing for the RSC 
has a good performance. Likely, from Fig. 5b, the errors of (V7.6 — VZ) and (igqret — igq) are about zero; in other 


words, the controller suggesting for the GSC has a qualified performance. 


When a measurement noise, 5% of rated values, is added to the measurement signals, i, w,, with the 
wind profile as Fig. 4a, the simulation results are shown in Fig. 6. This figure indicates that with the above 
noise measurement, the turbine using the recommended MPPT scheme still tracks its maximum points more exactly 
comparing with the turbine using the old MPPT scheme. However, comparing with the case of pure measurement as 
Fig. 4, the measurement noise causes a negative impact on the turbine performance but this impact is insignificant. 


Fig. 7 is the simulation results for a wind profile which varies rapidly as Fig. 7a. It is easy to see from 
Fig. 7b and Fig. 7c that the turbine using the suggested scheme has more qualified performance in both the power 
coefficient C, and the electrical energy in total comparing with the case using the conventional scheme. 
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Fig. 7. Simulation results in the case of rapidly-varying wind velocity: ( a) wind profile; ( b) power coefficient; ( c) 
electrical energy output 


6. CONCLUSIONS 

This paper suggested a new MPPT scheme for variable speed wind turbines using a DFIG. The suggested 
scheme allows the turbine tracking the maximum-power-point more effectively than the turbine using the MPPT- 
curve scheme. As can be seen from the simulation results, by using the suggested scheme, Cp almost keeps up 
Cymax» Wr Was approximate to wropt and the electric energy generated by the DFIG was higher comparing with that 
using the conventional scheme. Furthermore, by the suggested controllers, the rotor speed and current of the DFIG 
converged to their desired values. Thus, by the suggested control law, the wind turbine can achieve stable operations. 
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